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Abstract. - The apparent T 2 -temperature dependence of the London penetration depth in 
Sr2Ru04 is discussed on the basis of a multi-gap model with horizontal line nodes. The influence 
of the nodes in combination with nonlocal electromagnetic response leads to low-temperature 
behaviors as — T 2 \nT in the London limit and as T 2 in the Pippard limit. These behaviors 
appear only at very low temperature. On the other hand, the interplay of the superconductivity 
of different bands is responsible for the observed T 2 -like behavior over the wide temperature 
range. The experimental data can be fitted well with a set of material parameters as was used 
in the specific heat fitting. 



Since the discovery of the superconductivity in the layered perovskite compound S^RuCU 
, much effort has been devoted to the identification of its Cooper pairing symmetry 2 . An 
early theoretical prediction of the spin-triplet pairing 3 is supported by the NMR Knight 
shift measurement, which shows no change in the spin susceptibility on passing through the 
superconducting transition at T c (=1.5K) Moreover, /xSR measurements reveal the ap- 
pearance of an intrinsic magnetism with the onset of superconductivity, which was attributed 
to pairing with broken time reversal symmetry 0. These two experiments consistently iden- 
tify the superconducting phase as the chiral p-wave state di(fe) = zAi(T)(k x + ik y ), an 
analog to the A-phase of the superfluid 3 He. Here a representative form of the gap function, 
A a ^(fc) = (i<j v cr) a j3 ■ d(k), is given in the d-vector representation Such a gap function 

arises in a weak-coupling theory for a two-dimensional Fermi liquid with spin-orbit coupling 
013 El- This pairing state suggests a nodeless quasiparticle gap, Ai(T). 

As sample quality improved, however, power-law temperature dependence was observed 
in various quantities, such as specific heat UJ, NQR relaxation rate ^Oh thermal conductivity 
[111 I12j and ultrasonic attenuation 13.. These findings seemed to indicate the unexpected 
existence of line nodes in the gap and motivated the investigation of models with vertical 
(pseudo) line nodes ^JEDEj- Heat transport measurements, however, gave no evidence 
of strong anisotropy of the gap in the basal plane, essentially excluding vertical line nodes 
jlll 112) . If at all, nodes would have to be horizontal, i.e., parallel to the k x -k y -p\ane Still 
such a gap structure alone seemed insufficient to account for all the thermodynamic data. 

Motivated by these experimental situations and the fact that S^RuCU has three conduc- 
tion bands, called a, (3 and 7, Zhitomirsky and Rice (ZR) proposed a multi-gap scenario, in 
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which the primary nodeless gap, d\(k), is developed in the single "active" band (7-sheet), 
while the secondary gaps with horizontal line nodes of the form, 



appear in the other ("passive") bands (a, /3-sheets) through an interband Cooper pair scat- 
tering ^2]- Then this scenario was used to give a very good fit in the specific heat over the 
whole temperature range. The picture of an orbital dependent superconductivity (ODS) |18| 
provides a consistent view in many other aspects. Identifying the 7-band as the dominant 
band for superconductivity makes sense in view of the fact that ferromagnetic spin fluctua- 
tions, beneficial to spin-triplet pairing, are slightly enhanced in this band ^5JGI3]> while the 
other two bands show very strong spin correlations for the finite incommensurate wave vector 
[2HE21 In addition, the chiral p-wave state is stabilized by spin-orbit coupling relative to all 
other triplet states, if the 7-band is dominant |23) . 

On general grounds, however, all scenarios involving line nodes are hampered by the fact 
that line nodes can appear in a spin triplet pairing state only accidentally, if spin-orbit coupling 
is present. Blount has shown that, in general, only point nodes appear by symmetry reasons 
|24| . In order to generate rigorous line nodes the ZR scenario requires exclusively the interlayer 
pairing for the a-/?-bands. Any inplane pairing component would leave point nodes only, in a 
rigorous sense. Nevertheless, we will assume the presence of horizontal line nodes for the time 
being, and later discuss consequences of the presence of the inplane pairing components. 

In contrast to other thermodynamic quantities, the apparent T 2 -dependence over wide 
temperature range below 0.8 K observed in the London penetration depth, A(T), for fields 
along the z-axis, |25| seems incompatible with line nodes, since a simple-minded theory would 
lead a T-linear behavior. Bonalde et al. ascribed T 2 -behavior in their data to nonlocal 
effects, as was discussed by Kosztin and Leggett (KL) in the context of the high-temperature 
superconductors |25l |2E] . The KL-theory is based on the fact that for the large limit of the 
Ginzburg-Landau (GL) parameter k = A(0)/£o (London limit), a small portion of the Fermi 
surface in the vicinity of the nodes (~ k _1 ) should be treated in the nonlocal limit. This 
would alter the power-law from the linear to a T 2 -behavior. While this discussion is applicable 
rigorously for very large K and in the very low-temperature regime only, Bonalde et al. argued 
that the relatively small k (~ 2-3) of Si'2Ru04 could lead to a wider temperature range, even 
up to T* ~ A(0)k _1 ~ 0.6T C , below which the KL renormalization appears. Although there 
were several further theoretical explanations for the peculiar power-law behavior |23 |2Hj , 
a clear understanding has not been achieved yet. 

In this Letter, we show that the rigorous asymptotic power-law behaviors in A(T) due to 
nonlocal effects are restricted to very low temperatures in both the London (k ^> 1) and the 
Pippard (n <C I) limits. Instead, based on the ODS model the overall T 2 -dependence and 
beyond can be fitted well and consistently with the specific heat data. Thereby the non-locality 
of the electromagnetic response is essential for the fitting of the experimental data and the 
resultant Ginzburg-Landau (GL) parameter k = 2.3 is good agreement with the experimental 
value. 

We restrict our discussion to the clean limit and choose the sample geometry so that the 
x-z-plane surface is exposed to magnetic field parallel to the z-axis, so that both the vector 
potential A and the screening currents j are oriented parallel to the x-axis. The penetration 
depth is obtained from the electromagnetic response kernel, K(q,T) = —(4ir/c)j(q)/A{q), as 



d 2 (fe) = zA 2 {T)(k x + ik y ) cos(k z c/2), 



(1) 




(2) 



assuming a specularly reflecting surface. 
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We only consider intraband pairing so that in the lowest order each band separately con- 
tributes to the response kernel, 

*(«,T) = iZ;Cjtfj(«>r). (3) 



Here the zero-temperature penetration depth in the London limit has been defined as Ao = 
yjrac? /4nne 2 with n/m — ^2i( n l/ m l)i where ni and mi are the density and the effective mass 
of electrons in the band I. The ratio, £; = (ni/mi)/(n/m), gives a measure for the contribution 
of each band to the current density at T = 0. The dimensionless response kernel has the form 



K l (q,T) = 2TrT V V- / du - v = — , (4) 



2 r 1 , Vi - M^lrf/ 



n— — oo fc^ 







where w„ = (2n + l)7rT is the fermionic Matsubara frequency. Here, we have assumed cylin- 
drical Fermi surfaces for simplicity and used (d\di\ 2 / dk) • q = for the present geometry. 
The nonlocality appears via Qi = qv^i/2 = (jrq/2noi)Ai(0) where the dimensionless variables 
q = q\o and Kqi = Ao/£o; with the coherence length for each bands, £o; — vfi/kA.i(0), were 
introduced. For nodeless gaps or gaps with horizontal line-nodes, the angular integral /i in 
eq. (@J yields 



2 



where g(z) = (vl + z~ 2 — l)/z. At zero temperature, replacing 27rT^ n by /_ du>, we obtain 

^M)=5>(p|), (6) 

where Qo(z) = (2/z) tan _1 (z) — ln(l + z 2 )/ z 2 . 

We consider now the contributions from the bands with horizontal line nodes, i.e., the 
passive bands in the ZR scenario (I — 2), since their quasiparticles would dominate the 
low-temperature behavior, while those in the totally gapped dominant band are inactive. 
Analogous to the KL theory, we evaluate the average over k z for an approximate gap form in 
the vicinity of the nodes. It leads to 

5Ki(q,T) = Kt(q,T) - tf,(g,0) = 5Ki{Q,T)F{q/t), (7) 

where 6Ki(0,T) = — 21n2(T/A2(0)) is the usual T-linear contribution to the kernel and 
t = T/Tl with Tl = Q2/9 = (tt/2ko2)A2(0). For horizontal line nodes the universal function 
F(z) has a form different from the KL-expression: 

2^) = !-^ ^/(*)(l-f) 2 , (8) 

and f(x) = l/(e x + 1). Note that F(z) has the asymptotic form, F(z) sa 7r 2 /6zln2 for z ^s> 1 
and 1 — z/61n2 for z <C 1. Since K 2 {q, 0) = 1 — trq/Qn^ ~ 1 in the London limit (^02 3> 1) 
(eq. ©), the deviation of A(T) from its zero-temperature value, A(0) (= Ao), is given by 

AA(T) _ A(T) - A(0) _ 2 [°° d JK 2 (q,T) (g) 



Ao A 7T J (q 



4 
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An appropriate estimate of eq. © gives the following asymptotic expression for T < 0.05Xl, 

AA(T) it 2 [T\ 2 { T\ 

1 ' '1.0 + In— . (10) 



Ao 3k 2 \TlJ \ T L 

In the Pippard limit (kq2 <C 1) we use g(z) = z~ 2 in eq. (JjJJ. Then, the Matsubara sum 
can be evaluated straightforwardly 

k z 

From eq. J5J we obtain 



A(T) 
A(0) 



^ tanh' |d2 



7T \ - |d; 

2 ^ A 2 (0) """" V 2T 



-1/3 

(12) 



where A(0)/Ao = (27t/3a/3)(7t/4ko2) 1 ^ 3 . At low temperature the main contribution in eq. i|ll|) 
comes from k z ~ ±7r/c. Evaluating the fc z -average in eq. Ijlljl or using the large-g limit of 
eqs. © and JJJ, the low-temperature expression is given by 



AA(T) 1 / T \ 2 83 / T x 



A(0) 18 VTpy 3240 \T P 



(13) 



where Tp = A2(0)/7r. The deviation from the Independence becomes visible for T larger 
than 0.3T P . Note that the prefactor of the T 2 ter m in both limits is determined by the inverse 
of the characteristic temperatures, i.e. the smaller A 2 (0) gives the faster increase of A(T) 
close to T = 0. 

In conclusion, the power-law behaviors due to nonlocal effect are restricted to temperature 
below T*. For any k, T* has a value between both limits, i.e., 0.05Tl < T* < 0.3T P , which 
is much lower than A2(0). Since T* < 0.17T C even with the maximal gap, A 2 (0) = Abcs 
(= 1.76T C ), it already excludes a reasonable fitting of the T 2 -behavior with a single-gap model 
including horizontal line nodes for S^RuO^ 

We now turn to the multi-gap scenario for Si^RuCU 17 , in which the temperature de- 
pendence of the gap magnitudes are given by an effective two-band model; A a = Ap = 
A 2 (T) and A 7 = A\(T). In the gap equations, the pairing interactions are defined as 
Vxi = — gik ■ k 1 , V22 = — 2g 2 (fc • k') cos(k z c/2) cos(k' z c/2) for intraband coupling and V\2 = 
— \/2g^,{k ■ k') cos(k' z c/2) for the interband Cooper pair scattering. We estimate various other 
material parameters from the de Haas- van Alphen measurement |5]: the Fermi velocities, 
vp a r~j v-pp = 0.68t>F 7 , the (Fermi surface averaged) band masses over their carrier den- 
sities, Cq : C/3 : C7 = 0.20 : 0.44 : 0.36 and the densities of states at the Fermi level, 
^Voi : Ao 2 = 0.57 : 0.43. We fix the parameters as g\ — 0.4 (in units of the inverse total 
density of states), g 2 /.gi = 0.85 and ^02/^01 = (A 2 (0)/Ai(0))/0.68. This leaves two free 
parameters for fitting: gz/gi and Koi- 

Figure ^ shows that the temperature dependence of the gaps A/(T) for different interband 
couplings gz/gi- In the absence of interband coupling (93 = 0) a second phase transition 
occurs for A 2 at temperature T/T c ~ 0.2 yielding a gap of order A 2 (0)/Abcs ~ 0.2. The 
stronger the interband coupling the less visible the feature of this transition becomes. For 
strong interband coupling all gaps are tied together and act like a superconductor with a single 
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Fig. 1 - The temperature dependence of the gap functions scaled by the BCS gap at T — 0. The 
interband couplings are varied, while the other parameters are fixed as gi = 0.4 and gij g\ = 0.85. 

order parameter. The ZR fitting of the specific heat suggests that we are in an intermediate 
regime of interband coupling, (73/51 = 0.07. 

In Fig. [21 we compare the calculated penetration depth for various values of (73/51 and 
K01 = 2.0 with the experimental data by Bonalde et al. j^S] with £01 = 700A 29 . For the 
weak interband coupling, the penetration depth increases rapidly at low temperature below 
T/T c ~ 0.2, which reflects both the smallness of the passive-band gap, A2(0), and the rapid 
decrease of Aa(T) at T/T c ~ 0.2. For the strong interband coupling, the T 2 -increase due to the 
passive bands is easily overshadowed by the active-band contribution, showing a monotonous 




Fig. 2 - The temperature dependence of the penetration depth for the different ga/gi. The other 
parameters are fixed as gi = 0.4, gilg\ = 0.85 and «oi = 2.0. A moderate interband coupling fits 
well the observed temperature dependence. 
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increase with concave shape. Intermediate interband couplings give a temperature dependence 
with a weak convexo-concave shape by subtly combining the contributions from both bands. 
The experimental data indeed exhibit this type of behavior. The moderate interband coupling, 
93/91 — 0.07, which is exactly the same parameter as ZR to fit the specific heat measurement, 
also reproduces the observed temperature dependence over the wide temperature range. 




Fig. 3 - The koi dependence for gi = 0.4, (72/31 = 0.85 and gzjgx = 0.07. The resultant Ginzburg- 
Landau parameters, ki = A(0)/£oi are shown in the parenthesis. The thin solid line is for the nodeless 
gap in passive bands (see text in detail). 



Next we fix the interband coupling constant to the value optimized by ZR gs/gi = 0.07 
and use kqi as the only free fitting parameter. Figure shows curves for various values of 
kq\, where the resultant GL parameter for the active band, K\ = A(0)/£oi> is shown in the 
parenthesis. The slopes of the whole curve monotonously decrease starting from the London 
limit and approaching the Pippard limit. Irrespective of a measure of non-locality, all curves 
show the convexo-concave shape. The best fitting for the experimental data is given by 
Kqi — 2.0 with Ki = 2.3, which is slightly smaller than the generally cited value k ss 2.6 
|29j. We also estimate k%/ki = £oi/£o2 ~ 0.66. The small discrepancy between measured 
k and ours is very plausible, in view of the fact that k = 2.6 results from the measurement 
of the London penetration depth by means of /uSR in the mixed state. The order parameter 
in the passive bands with a larger coherence length is more easily diminished by a magnetic 
field than that of the active band. Consequently, the /^SR-measurement observes a larger 
London penetration depth, where screening currents from the a-/3-bands are reduced. This 
explains the smaller value of the very-low-field penetration depth as measured by Bonalde et 
al. jSHlEHj- Our fit suggests a low-field value of A(0) « 1600A compared to /i,SR-derived 
value of A(0) ~ 1900A which is rather likely only due to the 7-band contribution |29| . 

So far we have assumed horizontal line nodes in the passive bands due to interlayer Cooper 
pairing (ZR scenario). As an opposite limit, we consider the nodeless gap, d-i — z/S.2{T)(k x + 
iky), corresponding to the case of the intralayer pairing only. Trying the same fitting procedure 
we find only a minor difference at very low temperatures as shown in Fig. (thin solid line). 
While the difference is small in A(T), the removal of the nodes would be more easily visible 
in the specific heat data. 
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We have shown that the nonlocal response for the horizontal line node yields — T 2 lnT 
behavior in the London limit, while T 2 power law in the Pippard limit yet at quite low tem- 
perature. On the other hand, the ODS model in combination with nonlocal response naturally 
explains the apparent T 2 -dependence in a quantitatively consistent way with the specific heat 
fitting. Concerning the question of nodes, it seems that the London penetration depth data 
are less affected by nodes than the specific heat data. Eventually, a more anisotropic but 
nodeless gap shape would likely be sufficient for good fits in any case. In view of the compli- 
cated quasiparticle spectrum in the multi-gap compound Si^RuO,!, it is difficult from present 
experimental data to come to conclusive answer concerning the detailed gap shape. 
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